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ABSTRACT 



A reflex arc nitrogen plasma in a longitudinal magnetic field has 
been investigated by means of cylindrical Langmuir probes. Magnetic 
fields varying from 2400 gauss to 6600 gauss were studied. A radial 
profile of electron temperature and space potential was obtained as an 
aid to understanding the inner and outer rotational fluctuations previously 
noted to be characteristic of the reflex arc column. The effect of crossed 
magnetic and electric fields resulting in an E x B particle drift has been 
proposed as the primary mechanism of the outer rotational fluctuation, 
while the inner rotation mechanism has been proposed to be the j x B 
force of the magnetohydrostatic equation of a plasma in equilibrium. 
Diagnostic measurements of a theta pinch device on the plasma column 
were not successful due to circuit pick-up of electric signals due to 



the theta pinch current pulse. 
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I. INTRODUCTION 



This work is an examination of the basic parameters of a reflex 
arc, magnetically confined, nitrogen plasma column produced in the 
plasma facility of the United States Naval Postgraduate School. Of 
particular interest were the turbulent nature of the plasma column and 
the modification of previously-measured parameters as a result of 
installation of a vacuum column constriction and a larger diffusion pump 
on the plasma facility. In addition, the effect of the newly installed 
theta pinch coil was to be investigated. Langmuir probe studies were 
conducted to determine the characteristics of the plasma fluctuations 
and to furnish parameters to aid in establishing a theoretical basis for 
the oscillations observed. 

Langmuir probes and photomultiplier tubes were then to be used 
for plasma analysis during the application of the theta pinch and the 
expected theta pinch instigated instability period. The theta pinch 
device achieved operational status late in the period allotted for research 
and only a cursory look at the theta pinch effects was achieved. 



system, shown in Figure 1, consists of a four-inch diameter pyrex glass 
pipe connected to a high vacuum system capable of maintaining the 



II. THE PLASMA FACILITY 



The plasma facility of the Naval Postgraduate School has been 
described in detail by Streit and Olsen [l] and Andrews The 




The basic 
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Figure 1. Top View of Plasma Facility 






neutral gas pressure in the vicinity of 10 mm Hg while the system is in 
an operational state. The pyrex column is enclosed by six main magnet 
coils and a mirror magnet . Plasma -genera ting electrodes, a hollow 
cathode and annular anode, are located at the north end of the system 
and are contained in a complementary magnet system. The main magnet 
system is designed to produce a variable field up to 10,000 gauss which 
is homogeneous to within 2.5% along the axis of the plasma column. 

The cathode-anode circuit is shown in Figure 2. A negative 
pressure gradient aids the flow of the nitrogen through the cathode toward 
the vacuum chamber. Electrons thermionically emitted from the inside 
wall of the cathode collide with and ionize a portion of the nitrogen 
flow. The pressure gradient as well as the potential gradient between the 
cathode and anode draw the electrons that have not returned to the 
cathode wall out of the cathode and toward the anode. A highly visible 
beam of plasma can be seen emerging from the cathode and passing 
through the hole in the anode. Once through the anode the nitrogen 
column is in the reflex arc configuration. The floating anode is insulated 
and serves only as an end point for the plasma column. In this experi- 
ment a 60 ampere cathode current was maintained. This current 
represents a cathode-anode current as the reflex arc column has no axial 
net current. Argon was used to start the arc and to establish the initial 
plasma beam. As the argon was turned off, nitrogen was fed into the gas 
supply. This operation is readily observable in the plasma column as 
a slow shift in color from the blue of argon to the pink of nitrogen. 
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Figure 2. Cathode - Anode Arrangement 




The recently completed theta pinch device has resulted in the 
installation of a 2-inch diameter constriction in the 4-inch diameter 
pyrex pipe. The theta pinch device consists of a 6-inch long single 
turn copper ring which fits around the 2 -inch diameter pyrex glass 
constriction. The copper ring is connected in series with a 0.75 
microfarad capacitor bank. The configuration of the theta pinch device 
is illustrated in Figure 3. 

III. PLASMA FLUCTUATION STUDY 

A. BACKGROUND 

The plasma fluctuation studies were initiated with two objectives. 
The validity of the data obtained by Andrews LQ was to be re-examined 
in relation to the changes in the configuration of the plasma facility, 
and further studies were to be conducted into the characteristics of the 
plasma in order to interpret the initiation and propagation of the observed 
rotational oscillations. 

B. STUDY TECHNIQUES 

The techniques employed to confirm Andrews' findings used 
Langmuir probes, as opposed to Andrews’ use of photomultiplier tubes . 
Figure 4 shows the Langmuir probe installation configuration. The probes 
were situated in the ports between main magnets #5 and #6. The probes 
consisted of .02 inch diameter tungsten wire, shielded over part of their 
length by a stainless steel tube. The tungsten wire directly in the 
plasma column was surrounded by alumina, which also served as an 
Insulator Inside the stainless steel tube. The projecting probe length 
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Figure 3. Theta Pinch Device 
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Figure 4. Arrangement of Langmuir Probes Used to Confirm Rotational 
Direction of Oscillation 
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of bare tungsten wire was 0.117 inches. A number of different circuits 
were necessary for the conduct of the investigation and will be discussed 
in the following section. 

C. CONDUCT AND RESULTS 
1 . Introduction 

Previous investigations of the plasma produced in the facility 
in the reflex arc configuration have revealed three dominant periodic 
fluctuations. Andrews, through photomultiplier analysis , found what 
appeared to be two rotational fluctuations, termed inner and outer, with 
frequencies ranging from 15 kHz to 400 kHz, and a 60 Hz-360 Hz fluctu- 
ation which did not appear to be associated with any particular plasma 
motion. The rotational oscillations were described as having a counter- 
clockwise rotation when viewed from the cathode end of the plasma 
apparatus. It was the purpose of this investigation to confirm the 
presence of these fluctuations in view of the recently modified plasma 
facility parameters, and to obtain additional characteristics of the moving 
plasma column. 

2 . 60 Hz -360 Hz Oscillation 

The first approach to the 60 Hz - 360 Hz oscillation was to 
confirm the presence of the fluctuation while the plasma facility was in 
normal operation. The technique employed in this study was to allow 
the plasma column to reflect its fluctuation through a probe connected 
to ground through a 500 ohm resistance. Fluctuation in the plasma 
column was reflected as voltage variation across the 500 ohm resistor. 
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The arrangement shown in Figure 5, with the variable electronic filter 
adjusted to 400 Hz low pass, resulted in a 60 Hz sinusoidal oscilloscope 
trace with a 360 Hz signal riding upon the 60 Hz signal. The probe in 
the plasma and the dummy probe outside the pyrex tube displayed the 
trace in phase with one another. At this point, the dummy probe was 
moved from port to port toward the cathode end of the column with no 
phase shift being noted. The absence of a noticeable phase shift 
would show that the plasma column underwent the fluctuation at the 
same time along its entire length. It would have been difficult to detect 
wave propagation down the column if velocities were greater than 100 
meters per second. 

Two sources were considered likely candidates for the origin of 
this fluctuation; the cathode-anode arrangement and the magnetic field 
system, both of which are dependent on the public power system. 

The second phase of the investigation was conducted to determine 
if the magnetic field system was responsible for at least a portion of 
the 60 Hz - 360 Hz fluctuation. The plasma facility was brought to full 
operating conditions, short of establishing the gas flow and plasma 
column, with a magnetic field of 2400 gauss. Figure 6 is the resulting 
dummy probe oscilloscope display using the circuit shown in Figure 5. 

The 60 Hz - 360 Hz oscillation was clearly visible on the oscilloscope 
display tube both for the dummy probe and the probe in the vacuum 
column. The large undulating signal of Figure 6 is a 60 Hz wave, while 
the smaller signal riding on the 60 Hz signal is a 360 Hz signal. 
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Figure 5. Arrangement used to investigate 60 Hz - 360 Hz oscillation 




Figure 6 . 



Dummy probe oscilloscope display with no plasma and 
2400 Gauss magnetic field 
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The earlier photomultiplier tube studies conducted by Andrews 
show the presence of the 60 Hz - 360 Hz fluctuation in the light 
emission from the plasma column. This means that the fluctuation 
must be dealt with not only as an effect on the probes, which could be 
greatly reduced by using two probes with a differential amplifier, but 
also as an effect on the stability and characteristics of the plasma 
column . 

Calculations of plasma facility circuit parameters have also 
implicated the cathode-anode complex in the 60 Hz -360 Hz signal 
propagation, with the anode voltage variation being approximately 6% 
of the desired voltage setting. 

The origin of this signal, presumably the public power system 
acting through the cathode-anode voltage and the magnetic field system, 
could possibly be better controlled through filtering, but for the present 
this does not appear to be a practical possibility, primarily because of 
financial considerations. As the 60 Hz - 360 Hz signal cannot now be 
eliminated, the present study will involve working with or through this 
signal. 

3 . Radio-Frequency Inner and Outer Oscillations 

After initially determining that the previously observed inner 
and outer rotational oscillations could be registered on the Langmuir 
probes using the circuit shown in Figure 5, a study was made to confirm 
the direction of these rotations. The spaced probes as arranged in 
Figure 4 were used in conjunction with the circuit illustrated in Figure 7 
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Figure 7. Probe Circuit I sod to Determine iroction of ^^otation 
of Plasma Rotational Fluctuation 




Probe Fist 



Probe Upper 



Figure 8. Oscilloscope Display Reflecting Phase Difference of 
Rotational Fluctuation at Probe East and P’-obe Upper 
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to determine the direction of rotation. The periodic fluctuation of the 
plasma was observed as a voltage variation across the 500 ohm resist- 
ance in the probe circuit. The probes were positioned 3 cm from the 
center of the column and the phase relations of various probe combina- 
tions were examined. Figure 8, with probe east on upper beam and 
probe upper on lower beam, is an example of the dual beam oscillo- 
scope displays obtained. Results of the investigation agree with the 
previous photomultiplier tube studies of Andrews . The rotational 
direction was found to be counterclockwise when viewed from the 
cathode end of the plasma column. Similar observations of the inner 
rotation at a radial distance of 0.5 cm also show a counterclockwise 
rotation when viewed from the same position. 

Using the circuit of Figure 7, measurements were made at magnetic 
fields of 2400 gauss to 6600 gauss with the objective of establishing 
the plasma rotational frequency as a function of radial position. Once 
more the periodic plasma fluctuations were registered as voltage fluctu- 
ations across the 500 ohm resistance in the probe circuit. The results 
for the outer mode are shown in Figures 9 and 10 and for the inner mode 
in Figure 1 1 . 

The outer or slow rotation demonstrates two macroscopic properties 
of interest, other than the actual rotational velocities. As a general 
rule, for any given radial position within the plasma column, the fre- 
quency of rotation increased as the magnetic field increased. At the 
highest magnetic field examined, radial positions 1.0 cm and 2.5 cm 
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Figure 9. Outer Fluctuation Frequency as a Function of Radial Position and 
Magnetic Field 
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Figure 11. Magnetic Field versus Rotation Frequency for Inner Fluctuation 



violated this rule. A greater rotational velocity was found at the 1.0 cm 
radial position for all except the highest magnetic field investigated. 

These results are elements in the explanation of the propagation of the 
fluctuation and will be explored more fully in section III, paragraph D. 

The inner or fast oscillation was not observed at a radial position 
greater than 0.5 cm from column center. Figure 11 shows a nearly linear 
relationship between frequency and magnetic field over the range measured, 
with a slight deviation at 6600 gauss. 

Though the terms inner and outer have been applied to the rotations, 
Langmuir probes registered both phenomena at the 0.5 cm radial position, 
thus making it difficult to establish separate zones for the oscillations, 
if such zones exist. Movement of the probe into the column center 
resulted in increased random voltage fluctuations in the probe circuit, 
but outer oscillations were still noted. The outer and inner oscillations 
overlap within the center of the plasma column according to probe 
measurements. The 0.117 inch long probe may have been extending 
through two different oscillation zones of the plasma column. 

4 . Plasma Profile Studies 

An investigation of the plasma column in profile was under- 
taken to form a foundation for possible mathematical treatment. of the 
plasma fluctuations. Studies were conducted using a time sampling 
system similar to that used by Hart [3] , the objective being to measure 
the plasma parameters with a Langmuir probe during specific phases of 
the plasma oscillation. The circuit shown in Figure 12 was used to 
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Figure 12. Langmuir Probe Circuit for Plasma Column Time Sampling 





conduct this portion of the investigation. A -75 to +75 volt, 10 Hz 
sawtooth wave was impressed on the circuit as shown. Prior to the 
initiation of the plasma column, the resistance was adjusted so that 
a horizontal oscilloscope trace was obtained. Switch was closed 
only for testing the system, with resistance serving as the test 
probe. The dummy probe was placed adjacent to the vacuum column and 
test probe. Signals from the test probe and the dummy probe were passed 
into the differential comparator in the type 555 Tektronix dual beam 
oscilloscope where stray signals were removed and which passed on to 
the oscilloscope display tube the voltage differences between the test 
probe and the dummy probe. The x-axis of the oscilloscope trace is a 
direct display of the sweep voltage, -75 volts to +75 volts. The 
oscilloscope trace y-axis deflection amplitude is due to the voltage 
difference between the test probe and dummy probe. As the circuit 
diagram illustrates, the -75 volt to +75 volt sawtooth wave is impressed 
on both the test probe and dummy probe circuits at the same instant. 

As the dummy probe circuit has no path to ground except via the saw- 
tooth wave circuit and the differential amplifier, the dummy probe will 
assume the instantaneous potential of the sawtooth wave and no current 
will flow in the circuit. The test probe has a path to ground via the 
plasma and a current will flow through the 1780 ohm resistor. The 
vertical amplitude of the oscilloscope display represents the voltage 
difference between the dummy probe and the test probe due to the 
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potential drop IR in the probe resistor and can in turn be related to the 
current in the test probe circuit through the formula: 



E = 
R = 

^Test Probe 



I R 

1780 ohms 

, 

Vertical 
1780 ohms 



amperes 



The voltage drop across the 1780 ohm resistance is so slight that 
for measurement purposes the instantaneous sawtooth wave voltage may 
be assumed to be the test probe voltage along the x-axis of the oscillo- 
scope trace. The oscilloscope display is the probe characteristic curve 
of Langmuir's theory. 

F. F. Chen BD offers a comprehensive review of the theory of 

1 * \ 

the characteristic curve as shown in Figure 13. The following character 
istics apply to the labels used in Figure 13; 
i = Current to Probe 

Vp = Probe Voltage 

V = Space Potential (Probe same potential as plasma) 
s 

A = Region of Saturation Electron Current 

B = Transition Region or Retarding Field Region 

C = Region of Saturation Ion Current 

V^ = Floating Potential (Probe sufficiently negative to 
repel all electrons except a flux equal to the flux 
of ions) 
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Figure 13. Current-Voltage Characteristics of a Langmuir Probe 




Figure 14. Linear Extension of Ion Current Zone 
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The voltage applied to the probe is plotted.along the x-axis and 



the current flowing through the probe circuit is plotted along the y-axis. 
If the voltage applied to the probe is increased, three different parts 
of the probe characteristic will be traversed. In the region marked A in 
Figure 13, the current flowing to the probe is due to positive ions. 
Region B shows a rapid increase in electron current, which in region C 
reaches saturation. 

When a large negative voltage, relative to the plasma, is applied 
to the probe, the current is due to positive ions only, as the electrons 
cannot reach the probe. As the probe potential becomes less negative, 
the electrons in the tail of the Maxwellian distribution will begin to 
penetrate to the probe. 

A relation between the probe current (Ip), ion current (I.), and 

electron current (I ) can be written: 

e 

I = I -I- I. 
p e i 

In region B a Maxwellian distribution of electron energies is 

assumed to exist. If this is the case, then the fraction of the total 

electron current which reaches a probe at potential V is exp (-eV/kT ) 

e 

where T is the electron temperature, 
e 

The general expression for the probe current in region B becomes: 

I = I. + I exp(-eV/kT ) 
p i o e 

where I is the electron current in the absence of a retarding field, 
o 

Subtracting the ion current, the electron current is obtained. 

I - I. = I = I exp ( -e V / V T ) 
p i e o e 
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Further manipulation yields 



In I = In I - (e V / k T ) . 
e o e 

Thus, assuming the ion current decreases smoothly as the voltage 
becomes less negative, the electron current, calculated by subtracting 
the ion current (Region C) and its linear extensions from the measured 
probe current, can be plotted on a semilogarithmic scale versus the 
probe potential, as illustrated in Figures 14 and 15. The slope S on 
Figure 15 is equal to the value (e / k T^ ) . Knowing the slope of the 
graph, the electron temperature can be determined. 

The present study is concerned with the characteristics of the 
rotating plasma. The above described procedure presupposes the plasma 
remaining somewhat stable within the vicinity of the probe for at least 
the time of the voltage sweep. Two measuring methods seem to offer 
possibilities of overcoming the obstacle of a time dependent plasma. 

The first method would involve an extremely rapid voltage sweep of the 
probe. This method was not attempted for a number of reasons. Crude 
calculations reveal that to sweep during 1/20 of one cycle of a 30 kHz 
rotation would require a l.GAsec voltage sweep. A sweep of this speed 
would involve network impedance to a much greater extent than sweeps 
in the Hertz range. In addition, a study of pulsed probes by Bills, et al 
has shown that field-induced ion movement through a probe sheath 
may take as long as 1 microsecond. Most important, a sawtooth voltage 
sweep in the MHz range was not available with the voltage range desired. 
Initially, a 202 A function generator. was employed to furnish a triangular 
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Figure 15. Probe Characteristic on a Semilogarithmic Scale 




Figure 16. Schematic of Time Sampling Triggering Technique 
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wave, but the -18 volts to -t-18 volts maximum amplitude was found to be 
inadequate for the voltage characteristics of the plasma column. A 
definite ion current region was not visible under certain plasma condi- 
tions . This region was necessary for proper interpretation of the charac- 
teristic curve. The space potential of some regions of the plasma column 
was found to exceed 18 volts . 

The method adopted involved time sampling techniques similar to 
those employed by Hart, as mentioned earlier. The method employs the 
plasma oscillation as a trigger for a pulse which in turn modulates, the 
intensity of the oscilloscope beam, thus recording the characteristic 
curve only when the periodic plasma fluctuation is at a specific stage. 
The sampling probe of Figure 12 and the illustrated circuit were used for 
this purpose. As an aid to sampling, the trigger pulse could be delayed 
for a specific variable time so that various points of the oscillation 
cycle could be sampled. A schematic of the wave form seen by the time- 
sampling probe is shown in Figure 16. The radio-frequency oscillation 
superimposed on the 360 Hz fluctuation was used to obtain the triggering 
point. The point was selected at the minimum vertical displacement of 
the trace because of the better definition of the trace at that point. 

The probe sawtooth voltage sweep rate used was 10 Hz, which 
would result in 36 oscilloscope beam modulations per sawtooth sweep if 
the triggering took place once at the minimum vertical point of the 
360 Hz signal. The results of this technique offered a substantial 
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improvement in data resolution over the unmodulated probe characteristic 
curve. Figure 17 illustrates the effect of beam modulation, the lower 
two sweeps being unmodulated traces, and the upper display of dots 
being an unrelated time sampled sweep. 

Using the time sampling technique as illustrated in Figure 12, a 

radial analysis of the reflex arc plasma column was conducted from 

column center to d 5.0 cm radius, with fields of 2400, 3600, 4200, and 

4800 gauss. Anode voltage was in the order of 70 volts and the cathode 

was grounded. The cathode-anode current was maintained at 60 amperes 

Neutral gas pressure varied somewhat, with an attempt being made to 

obtain lowest possible background gas pressure at all times. The 

neutral gas pressure, as measured with an ion gauge, was generally in 

-4 

the vicinity of 10 mm Hg. 

The Langmuir probe characteristic curves were first analyzed for 
the space potential of the plasma column. The space potential versus 
delay values "as obtained from Polaroid photographs of the time sampled 
traces for a magnetic field of 2400 gauss are illustrated in Figure 18. 

At the radial distance of 0.5 cm, the time-sampling probe was unable to 
perform its function because of the presence of random fluctuations in 
the plasma at that point. As a result, the space potential for that 
position was obtained from a series of unmodulated traces. 

The periodicity of the space potential agrees with the earlier 
findings of periodic rotation of a portion of the plasma column. Too 
strict a value must not be attached to the apparent time of oscillation 
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H Figure 17. Z-axis Modulated 

(upper) and Unmodulated 
(middle and lower) Current 
versus Sweep Voltage 
Characteristic Curve 




-~6 

Delay (x (20) (10) seconds) 

Figure 18. Time-Sampled Plasma Space Potential as a Function of 
Radial Position 
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as reflected in Figure 18, as the readings noted were taken only at 
specific delay intervals, with the lines between points being added for 
clarity only. Missing points cannot readily be deciphered from the 
photographs and points for R = 1.5 cm and R = 2.0 cm have been deleted 
from Figure 18 for clarity. 

A compilation of data, similar to that illustrated above, was made 
for the reflex arc column in magnetic fields of 2400, 3600, 4200, and 
4800 gauss. The rotational oscillation study showed that different 
radial positions in the plasma column reflected differing oscillation 
frequencies. This situation would allow the averaging of space potential 
variations found for any particular radial position if it was found desirable 
to treat a particular radial position as having a specific space potential. 
The time sampled space potential measurements showed a time dependent 
space potential variation which would indicate a local azimuthal electric 
field and a variation of space charge density in the plasma column. 

A proposed mechanism illustrating the apparent situation is shown 
in Figure 19. It is difficult to compare this model with a macroscopic 
physical system as the rotating shells are not azimuthally symmetric with 
respect to charge. The periodicity of the phenomena at each radial 
position allows the meaning of the word shell to be stretched slightly 
to include the varied plasma parameters found at each radial position. 
Cylindrical symmetry is not implied. When consideration is given to 
the radial dependent azimuthal drift velocities, reference to an azimuthal 
electric field becomes meaningless, as in effect there would be a series 



38 



Density fluctuations 




Figure 19. Proposed Physical Mechanism of Rotational 
Fluctuation 



39 



of random, ever changing , radius -dependent azimuthal electric fields. 

A further hypothesis in explanation of the time dependent space 
potential variation would assign the establishment and propagation of 
the rotational oscillation to a space-charge-dependent effect, independ- 
ent of the effect of the sign of particle charge, with the electrons by 
virtue of their light mass leading the heavier ions which would have a 
larger frictional force with which to contend. The time dependent space 
potentials measured for each radial position were averaged to obtain an 
average space potential for the radial position. This allowed compari- 
son of the space potentials at various radial positions. 

Figure 20 is a graphical representation of the results of this 
comparison. The graph shows the presence of a radial electric field *E^. 
As the magnetic field- is increased, in the vicinity of the column 
center also appears to increase. Experimental results show that the 
column center remains in the vicinity of ground potential. This result 
would be expected, as the cathode remains at ground. The interior of 
the plasma column could be assumed to be fairly well shielded from the 
anode and would act as an extension of the cathode. 

The electron temperature distribution as obtained from the 
characteristic curves is shown in Figure 21. The point at 3600 gauss 
and 2.5 cm appears to be inconsistent when compared with the majority 
of the data. The electron temperatures indicated appear to be somewhat 
higher than the electron temperatures obtained by Cote using 
spectroscopic techniques. Cote's values being; 
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Figure 20. Space Potential versus Radial Distance as a Function of Magnetic Field 
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Figure 21. Electron Temperature versus Radial Position as a Function of Magnetic Field 
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These space-averaged electron temperatures are indicated by a bar 
graph adjacent to the space-resolved electron temperatures of Figure 21. 

Both the spectroscopic and Langmuir probe techniques show an 
electron temperature peak at the 4200 gauss magnetic field. The 
experimentally obtained electron temperature distribution indicates the 
presence of an electron temperature gradient in the plasma column, 
the high electron temperature being at column center. 

Gall and Oleson conducted a study of a steady state argon 
plasma in a longitudinal magnetic field by means of cylindrical 
Langmuir probes. The Naval Postgraduate School plasma facility was 
used in the same configuration as the present study, minus recent 
changes, with a hollow cathode discharge serving as the plasma source. 
Electron temperatures as a function of radial position were determined 
in magnetic fields ranging from 600 gauss to 9000 gauss, and are shown 
in Figure 22. The argon and nitrogen electron temperature profiles are 
quite similar, with the argon column core temperatures being higher 
than the nitrogen temperatures. At the 1.0 cm radial position and out- 
ward, the electron temperatures have similar values, within experimental 
error, for both argon and nitrogen. 
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Figure 22. Electron Temperature versus Radial Position as a Function of Magnetic Field 
for Argon 



D. INTERPRETATION OF RESULTS 



I . Space Potential Study 

The space potential study indicated the presence of a radial 
electric field directed toward the discharge centerline. As the rotational 
direction of the observed oscillation was in a right-handed direction 
with respect to the magnetic field, an*^x B effect is strongly implicated 
in the propagation of the rotational fluctuation. 

Representing the situation mathematically, the frequency of rota- 
tion can be described by the equation 

V 



Assuming the radial electric field observed to be perpendicular to the 
axial magnetic field, the classical solution to the motion of charged 
particles in such a crossed field situation yields a drift velocity perpen- 



-4 “4 

dicular to the crossed fields in the Ej-xBq direction with a drift velocity 

V , = E/B (2) 

d r o 

which is charge sign-independent. The drift velocity would assume an 
azimuthal form in the case under examination. In effect 

V = V (3) 

o d 

Combining equations (1), (2), and (3) above. 



f = 



2 TT r B 



(4) 



o 



Equation 4 indicates an inverse proportionality of frequency to magnetic 
field, opposite to the experimental results. Offsetting the inverse 



proportion effect of the magnetic field is the radial electric field, which 
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increases with magnetic field. The nature of the reflex arc is of 
interest at this point. The plasma column core apparently retains the 
potential of the cathode, while as radial distance in the plasma column 
increases the anode assumes the role of end plate for the plasma 
column. The plasma column becomes a portion of the circuit between 
cathode and anode. 

2 . Plasma Column Current Flow 

The space potential distribution would result in current 
flow, across magnetic field lines, from the maximum space potential 
region to column center. There would be no net current to the walls of 
the container as the insulated vacuum column walls could not continue 
to build up charge indefinitely. The proposed mechanism of current 
flow from the region of maximum space potential to column center is 
electron movement from column center to the region of large positive 
potential and hence to the anode. The cathode-anode current is 
considered to be the source of replenishment necessary for the steady 
state conditions achieved by the plasma column. No net current down 
the plasma column is noted. The higher electron temperatures at column 
center demonstrate high electron velocities at column center and add 
evidence in favor of a radial, cross-magnetic field current which the 
principles of thermal diffusion require. Exact calculation of this factor 
would involve knowledge of the distribution of the particle densities and' 
energies, which have not yet been determined. 
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A cross -field current suggests the possibility of a 



?xt 



force , 



which will also be in a right-handed direction with respect to the 
longitudinal magnetic field. This force is examined by Jancel , et al, 
in their discussion of the magnetohydrostatic equations. The 
magnetohydrostatic equations are obtained when the Maxwell equations 

V X B = 4 7T j and V • B = 0 

are added to the equations of plasma dynamics . 

The plasma equation of motion in equilibrium configuration 



becomes 



7 p =txt V u 



in which 



p = pressure 
B = magnetic induction 
P = mass density 

U = gravitational potential (negligible in 
magnetic field) 



The pressure term represents the force on the plasma. A state of constant 
oscillation is achieved. This implies that the force described above is 
opposed in some manner. Neutral particle pressure and velocity- 

dependent diffusion rates are proposed as mechanisms limiting the 

— ^ ^ 

particle motion resulting from the j x B force. 
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3 . 



Reflex Arc Oscillation Sources 



The E X B effect and the j x B force are proposed as the 
mechanisms by which the reflex arc plasma column reaches an oscil- 
latory steady state condition. This theory agrees in part with that of 
Datlov, et al, j^Q^ .who conducted a study of the characteristics of a 
reflex discharge (EIQ) in hydrogen in a magnetic field. The depend- 
ence of the origin of low-frequency rotational fluctuation on the 
parameters of the discharge was studied. The radial profile of the 
plasma potential was measured with a heated probe. The fluctuation 
frequency was found to be approximately proportional to the drift 

/B computed^ from the measured radial field 



velocity . v , = 
d 



(e 

V radial 



and magnetic induction. The fluctuation frequency was found to 
increase with magnetic field, but the growth of the frequency was 
explained by the fact that E^ increased somewhat more than linearly 
with the magnetic field B. 



In considering the results of Datlov with the results of this 



— ^ "4 

investigation, it is further proposed that the E x B effect mainly 
concerns the outer fluctuations and that the j x B force can be 
considered to be the primary effect involved with the inner oscillation. 
The linear dependence of the inner rotation frequency on magnetic 
field supports this proposal. 

The positive column offers an interesting contrast to the reflex 
arc column. The literature contains many references to the screw 
instability of a positive column, the theory of Kadomtsev C^Oj being 
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among the more widely noted treatments of the subject. Hoh and 
associates have contributed a number of important papers 

on the positive column. The key difference between the reflex arc and 
the positive column is the presence of a longitudinal electric field in 
the positive column. The instability mechanism is generally considered 
to be an off axis current filament drift which creates a charge separation, 
the ions assumed unaffected by the magnetic field lagging behind the 
rotation of the electron distribution. The space charge distribution in 
turn results in an x B drift as illustrated in Figure 23. This drift 
tends to force the plasma column toward the container walls. The axial 
electric field tends to turn the electron screw around the tube axis . 
Because electric neutrality is to be maintained, the ion screw must 



rotate with the electron screw. The role the j x B force plays in the 
phenomena is not yet clear. 

The mechanism for the reflex arc model, while employing an 
X B drift, does not appear to incorporate the Ey/-x B drift in the 
observed fluctuation. The differing drift velocities at different radii 
in the plasma column do not appear to support the establishment of 
an elec trie field. 

4 . Modification of Plasma Characteristics and Parameters 

The recent modification of the plasma facility created no 



apparent changes in the plasma characteristics. 
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Figure 23. A Cross-sectional View of Space Charge Distribution 
Responsible for Screw Instability of Positive Column 
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E. RECOMMENDATIONS 

1 . 6 0 Hz - 360 Hz Oscillation 

The feasibility of design of a feedback system to suppress 
the 60 Hz - 360 Hz fluctuations should be examined- The most probable 
employment of such a device would be to modulate the anode voltage 
supply so as to suppress the 60 Hz - 360 Hz oscillation imposed on the 
plasma column by both the cathode-anode complex and the magnet 
system. 

2 . Rotational Fluctuations 

The reflex arc plasma column should be examined for an 
onset and cut-off of the rotational oscillations and, if found, a study 
of plasma parameters in the vicinity of these magnetic fields should be 
conducted with the objective of establishing the conditions which lead 
to the onset of fluctuation and the continued propagation of the fluctua- 
tions . A study of plasma parameters using plasmas of other elements 
may prove beneficial in establishing the causes of rotational modes. 

The role of the anode in the propagation of the fluctuations should be 
investigated by varying the anode parameters. Possible areas, of investi- 
gation could include anode configuration and the effects of shielding 
the anode from the plasma column. 
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IV. THETA PINCH PERTURBATIONS 



A. BACKGROUND 

One of the devices that has been considered for possible use in 
a thermonuclear reactor has been the plasma pinch. The theta pinch 
device installed in the plasma facility is designed to compress a portion 
of the plasma column by a sudden increase of longitudinal magnetic 
field. This field increase is caused by a current pulse on the copper 
conducting ring around the plasma column. The column compression 
results in a small volume of high temperature, high density plasma. 

The pinch effect enjoys only a brief existence due to a rapid breakup 
once the pulsed magnetic field is removed. As noted by Rose and 
Clark DU , the imploding magnetic wall consists of lines of 
induction, the z-axis being the longitudinal axis of the vacuum column. 
The increasing B^ in turn induces a current i around the plasma 
column. Operational aspects of the device are illustrated in Figure 24. 

The installation of the theta pinch device was completed just 
prior to the completion of the period available for research. The device 
was tested with nitrogen gas in the plasma system. At pressures 
ranging from 0.5 microns to 80 microns, activation of the theta pinch 
device produced ionization in the nitrogen gas and propagation of an 
ionizing front down the pyrex column. These effects were visually 
observed . 
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Side View - Magnetic Induction Increases 




End View: A. Increasing Magnetic Field induced i to 
_^oi^ter rise in field. 

B. i X B force exerts pressure on column, 
compressing column. 



Figure 24. Theta Pinch Operation 
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REFLEX ARC NITROGEN PLASMA 



Attempts to observe the effects of the theta pinch operation on 
the reflex arc nitrogen plasma in a 1300 gauss magnetic field were at 
first unsuccessful. Photomultiplier tubes aligned on the plasma column 
and probes within the column both failed to register a noticeable change 
in the plasma characteristics. Upon each theta pinch discharge a great 
deal of electric noise was registered by both the electric probe and the 
photomultiplier tube circuitry. This effect could have hidden the theta 
pinch' effect on the plasma column, but if such effects were present, 
they would have been quite small compared to the normal plasma activity. 
Recent studies by plasma laboratory personnel have shown the presence 
of plasma perturbation. Further studies are in progress. 

C. RECOMMENDATIONS 

1 . Plasma Facility 

The theta pinch device should be modified to produce a greater 
pinch effect by the addition of a larger capacitor bank. Various plasma 
facility parameters, to include magnetic fields and neutral gas pressure, 
should be adjusted in the search for the effects of the theta pinch 
device, and an effort should be made to reduce plasma fluctuation. 

2 . Measuring Devices 

Probe and photomultiplier tube circuitry must be redesigned 
with shielding from the electric noise of the theta pinch device being of 
primary consideration. The use of light pipes should be investigated as 
an aid to removing the photomultiplier circuitry from the vicinity of the 
electric disturbances, 
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V. SUMMARY 



The use of cylindrical Langmuir probes in a reflex arc, magnetically 
confined nitrogen plasma confirmed the presence of fluctuation in the 
plasma column. The fluctuations were identified as being a 60 Hz - 360 Hz 
oscillation imposed by the power supply to the plasma facility, and two 
rotating fluctuations present in the outer and inner zones of the plasma 
column . 

Analysis of electric signals and the characteristic curves of the 
Langmuir probes had the following results: 

1. The 60 Hz - 360 Hz oscillation was found to be a character- 
istic of the plasma column as a whole with no wave propagation down 
the plasma column being noted. 

2. The inner and outer rotational fluctuations were shown to be 

in a right-handed direction with respect to the longitudinal magnetic field. 

3. The space potential profile of the plasma column appeared to 
resemble a gauss ian distribution with the column center at zero volts 
amplitude and the point of maximum space potential moving toward 
column center as the longitudinal magnetic field was increased. 

4. The electron temperature distribution profile also resembles 
a gauss ian distribution with the maximum temperature located at column 
center and the electron temperature decreasing with radius. 

The direction of rotation, space potential distribution, and electron 

i 

temperature profile lend credence to the hypothesis that an E x B drift 

^ 

effect and a j x B force are responsible for the rotational fluctuation 



propagation . 
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